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The formalism of non-standard four-fermion interactions provides a convenient, model- 
independent way of parameterizing a wide class of "new physics" scenarios. In this article, we 
study the performance of reactor and superbeam neutrino experiments in the presence of such non- 
standard interactions (NSI). Due to interference between the standard and non-standard amplitudes, 
sizeable effects are to be expected if the NSI parameters are close to their current upper limits. We 
derive approximate formulas for the relevant oscillation probabilities including NSI, and show how 
the leading effects can be understood intuitively even without any calculations. We will present a 
classification of all possible NSI according to their impact on reactor and superbeam experiments, 
and it will turn out that these experiments are highly complementary in terms of their sensitivity 
to the non-standard parameters. The second part of the paper is devoted to detailed numerical 
simulations, which will demonstrate how a standard oscillation fit of the mixing angle #13 may fail 
if experimental data is affected by NSI. We find that for some non-standard terms, reactor and su- 
perbeam experiments would yield seemingly conflicting results, while in other cases, they may agree 
well with each other, but the resulting value for #13 could be far from the true value. This offset 
may be so large that the true #13 is even ruled out erroneously. In the last section of the paper, we 
demonstrate that reactor and superbeam data can actually establish the presence of non-standard 
interactions. Throughout our discussion, we pay special attention to the impact of the complex 
phases, and of the near detectors. 
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1. INTRODUCTION 

At the dawn of the era beyond the standard model, a 
plethora of new theoretical models has been devised to 
resolve many of the experimental and theoretical short- 
comings of our current picture of elementary particles. 
However, in the context of future experiments, it is of- 
ten desirable to describe new physics in a more model- 
independent way. One possibility to achieve this is 
through effective four-fermion operators, so-called non- 
standard interactions (NSI), which arise naturally in 
the presence of heavy mediator fields. In this arti- 
cle, we shall focus in particular on NSI in the neu- 
trino sector, which have been discussed ongeneral phe- 
nomenological grounds in [T|, [H, 0, 0, M, H E , Hi , and 
in the context of specific models in [£J, [ToT [ill . [l2T ] . 
The importance of NSI for neutrino oscillation physics 
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has been pointed out in a pioneering work by Gross- 
man Qjl, and many _authors have studied their impact 
ITU [Till . \vi\ , atmospheric neutri- 
23p , conventional and upgraded 
26 



on solar neutrinos 

nos [n, M M Si 




30( 1 . neutrino fac- 
beta beams [11], 
cosmological relic neu- 
43j , neutrino-electron scatter- 
and neutrino- nucleus scattering [H, [Hi ■ Exist- 



neutrino beams 
tories @, [H, M 
supernova neutrinos 
trinos [Hi, e + e~ colliders 
ing [ 

ing experimental bounds are presented in [4 

Our main interest in this work will be on non-standard 
interactions in upcoming reactor and accelerator neutrino 
experiments. Although the main design goal for these ex- 
periments is the precision measurement of the standard 
oscillation parameters, the search for deviations from the 
standard framework is an equally interesting part of their 
physics program. Moreover, while in the race for the 
standard oscillation parameters, reactor and beam ex- 
periments are competing, we will show that their results 
will be highly complementary when one is interested in 
non-standard physics. 

In the numerical simulations which we are going t o 
present, we will focus on the experiments T2K [4a, 143] . 
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NO^A [50], Double Choo z I5ll [53 ] . and a hypothetical 
200 t reactor experiment [53|]. Of course, the analytical 
results apply also to other experiments such as T2HK |4§ ] 
and Daya Bay [H3 ]. 

We will first introduce our formalism in Sec. El and 
give a detailed discussion of the different possible Lorentz 
structures and their relevance for reactor and superbeam 
experiments. Although this discussion may seem rather 
technical, it will ultimately allow us to greatly simplify 
the problem and considerably reduce the number of pa- 
rameters. In Sec. [31 we will then present approximate 
expressions for the oscillation probabilities including NSI 
for the v e — ► v ei — ► z^ e , and — ► channels. We will 
also show in an intuitive way why certain NSI terms ap- 
pear in these expressions, and others do not. Sec. [His de- 
voted to a discussion of numerical simulation techniques, 
and of the specific experiments which wc have simulated. 
In Sec. 03 we show how the data from these experiments 
may be misinterpreted, if NSI are not taken into account 
in the fits. We will finally demonstrate in Sec. [6l that 
a combined analysis of reactor and superbeam data may 
allow for the actual discovery of a wide variety of non- 
standard interactions by goodness-of-fit arguments. Our 
conclusions will be presented in Sec. [7] 



2. THE FORMALISM OF NON-STANDARD 
INTERACTIONS 

2.1. The NSI Lagrangian 

It is well known that in the low energy regime, weak 
neutrino interactions can be described by effective four- 
fermion operators like 

C = ^ [P Q7 p (l - l 5 y a ] [7'7p(1 - 7 5 )/] , (1) 



and 



£msw = ^§ [v a Y{l ~ 7 5 H] [/7p(1 ~ 7 5 )/] , (2) 

where v a is the neutrino field of flavor a, £ a is the corre- 
sponding charged lepton field, and /, /' are the compo- 
nents of an arbitrary weak doublet. 

The low-energy fingerprint of many "new physics" sce- 
narios has a structure similar to Eqs. JT]) and (0), and 
the corresponding operators are called non-standard in- 
teractions. If we consider only lepton number conserving 
operators, the most general NSI Lagrangian reads 



£nsi = £>v±A + £s±p + £t, 



(3) 



where the different terms are classified according to their 
Lorentz structure in the following way: 



% J2 gS 4' f ' V±A l^ 1 - 7 5 )^] [/'7p(1 ± 7 5 )/] ■+ 
G| Yfflt*** [* a7 >(l - 7"H[/7p(1 ± 7 5 )/] + h. 



(4) 



C s±p = ^ E ' S±P [^(1 + 7 5 )^] [/"'(I ± 7 5 )/] , 
VZ fj' 

(5) 



C T = -4 E T \^ PT ^\ U'oprf] ■ (6) 



Here, Gf, is the Fermi constant, v and £ are the neu- 
trino and charged lepton fields, and the f's represent the 
interaction partners of the neutrinos. The dimcnsionlcss 
parameters e give the strength of the non-standard in- 
teractions relative to Gf, where an upper index s stands 
for NSI in the neutrino source or detector, while m de- 
notes non-standard matter effects, i.e. NSI affecting the 
propagation. In general, the s s can be arbitrary complex 
matrices, while the e m have to be hermitian. 

Note that we have required the neutrino fields to 
be purely left-handed, since processes involving right 
handed neutrinos would require either a neutrino helic- 
ity flip, or their amplitudes would have to contain at 
least two NSI terms (e.g. one to create the right-handed 
neutrino and one to absorb it), and would therefore be 
strongly suppressed. This constraint on the neutrino chi- 
rality in particular forbids vvjj terms in Cs±p and Ct- 

Before proceeding, let us give a simple estimate which 
relates the magnitude of the e parameters to the corre- 
sponding new physics scale Mnsi [H[ : If we assume the 
non-standard interactions to be mediated by some inter- 
mediate particles with a mass of order A/nsi, the effective 
vertices in Eqs. (J4j) ((6j) will be suppressed by 1/A/^ SI 



in the same way as the standard weak interactions are 



suppressed by 1/M-JjL Therefore we expect 



i£w 



(7) 



C 



V±A 



2.2. Relevance of the different NSI terms to 
reactor and superbeam experiments 

We see from Eqs. ((U) - © that the number of pos- 
sible NSI terms is very large. However, the number of 
parameters for our discussion of reactor and superbeam 
experiments can be greatly reduced by a few simple, but 
rather technical arguments. Many of these arguments 
are based on constraints coming from the requirement 
of interference between the standard and non-standard 
amplitudes. Of course, the total interaction rate will 
also contain pure NSI terms, for which these constraints 
do not apply; but they are suppressed by e 2 , and can 
therefore be assumed to be negligible compared to the 
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interference terms, which are linear in e. The following 
arguments are also summarized in Tab. |TJ 

1. The standard production and detection processes of 
reactor and superbeam neutrinos are, on the fun- 
damental level, decays of u quarks into d quarks, or 
vice-versa. Since interference of standard and non- 
standard amplitudes requires the external particles 
to be identical, only the e s ^^ terms with / = it, 
f = d will be relevant, and we will henceforth sim- 
ply omit the indices / and /'. 

2. For the non-standard matter effects, only coupling 
to electrons, up quarks, and down quarks is impor- 
tant. 

3. Non-standard couplings involving r leptons are ir- 
relevant since r production is impossible in reactor 
and beam sources, and is not considered as a de- 
tection process here, although it might in principle 
be possible for high energy superbeam neutrinos. 
Hence we take 



~s,V±A 



~s.S±P 

£ T0 



~s,T 







(8) 



For the same reason, processes involving muons can 
be neglected in reactor experiments, and processes 
involving electrons can be neglected in the super- 
beam source, since they constitute subdominant 
backgrounds even in the standard framework. 

4. In the couplings to muons, there is still room 
for non-(V — A)(V — A) contributions. For neu- 
trino production in pion decay, the effect of 
(S + P)(S ± P) type NSI is even enhanced by a 
factor of [5a, \5& 



rrif, m u + m d 



20, 



(9) 



and the importance of this enhancement for accel- 
erator neutrino experiments has been pointed out 
in [26| . However, there exist limits on the muon he- 
licity in pion decay [57], HI], which ensure that, in 
spite of the enhancement, (S + P)(S ± P) type NSI 
cannot affect the neutrino oscillation amplitude by 
more than a few per cent. 

5. Tensor interactions are impossible in pion decay 
since the decay operator must have a parity-odd 
component. 

6. In the detection processes involving muons, the 
(S + P)(S ± P) and TT terms are chirally sup- 
pressed by the smallness of m M compared to the 
typical superbeam energies of 0(1 GeV). As men- 
tioned above, the leading effect in the total event 
rate is given by the interference of the non-standard 
amplitude and the standard (V — A)(V — A) am- 
plitude. This interference can only occur if the ini- 
tial and final state particles have identical hclicitics, 



so for (S + P)(S± P) and TT type non-standard 
interactions, a mass-suppressed hclicity flip of the 
muon is required. We can also see the emergence 
of the suppression factor explicitly by considering 
the Dirac traces, which have to be evaluated when 
calculating the cross section. For example, in the 
case of (S + P)(S + P) NSI, the spin sum in the in- 
terference term of standard and non-standard am- 
plitudes is 

^Tr[ 7 "(l- 7 5 V / 2(l- 7 5 )^] 

spins 

•Tr[ 7p (l- 7 5 )^(l- 7 5 )dd] 
= Tr [7 P (l - 7 5 )(?V + «V)(1 - 7 5 )(^ + mu)] 

spins 

• Tr[ 7p (l - 7 5 )(^ + m M )(l - 7 5 )(^ + m d )] . 

(10) 

Similar equations can be derived for (S + P)(S — P) 
and TT interactions. Due to the orthogonality 
property of the chirality projection operators, a 
contribution proportional to remains of the first 
trace in Eq. (|10[) , and a contribution proportional 
to m u from the second. This leads to a suppres- 
sion factor of 0{m fJll m u /E 2 ). Low energy neutrinos 
(E < 1 GeV) interact with whole nucleons, rather 
than single quarks, therefore m u should be replaced 
by the much larger nucleon mass m n , so that in 
this case, the overall chiral suppression is only of 
0(m^/ E). At typical superbeam energies around 
1 GeV, we are in the transition regime between 
neutrino-nucleon interactions (quasielastic scatter- 
ing and resonance scattering) and neutrin o-q uark 
interactions (deep-inelastic scattering) (59T |60| . 

7. For (V — A)(V + A) interactions involving muons 
in the detector, chiral suppression occurs only for 
the hadronic interaction partners, and according to 
our above discussion, it is not very pronounced for 
them. Therefore, (V — A)(V + A) type interactions 
may in general be important for the cross sections, 
and modify their overall magnitude as well as their 
energy dependence. 

8. From measurements of the electron angular dis- 
tribution in nuclear (i decays, (S + P)(S ± P) 
and TT couplings to electrons are strongly con- 
strained [oil . wX . HH, l&il . |65| . Consequently, we take 



~s.S±P 
£ e0 



~s,T 
£ e& 







(11) 



for P 



9. There is still room for (V — A)(V + A) type terms 
involving electrons, because these terms differ from 
the standard model term only in the quark current, 
which cannot be directly measured. Limits exist 
only for the effective vector and axial-vector cou- 
plings of protons and neutrons [65[ , but due to the 
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Reactor source and detector (/ = u, f = d) 





£ a — e 


Source 
£a = H 




£ oc — 6 


Detector 

£a = H 


£ a =r 


V 


-A 




/ 


no fj, production 


no r production 


/ 


no n production 


no 


r production 


V 


+ A 




/ 


no fi production 


no t production 


/ 


no n production 


no 


t production 


s 


- P 


stronj 


; constraints 


no (i production 


no r production 


strong constraints 


no n production 


no 


r production 


s 


+ P 


stronj 


I constraints 


no /i production 


no t production 


strong constraints 


no n production 


no 


t production 


T 




stronj 


I constraints 


no fi production 


no r production 


strong constraints 


no n production 


no 


t production 



Superbeam source and detector (/ = u, /' = d) 





£a — 6 


Source 

£a = H 


£ a =T 


£a — £ 


Detector 
la = M 


£ a =T 


V -A 


no e production 


/ 


no t production 




/ 


/ 


no 


r detection 


V + A 


no e production 


/ 


no r production 


/ (mild supp.) 


/ (mild supp.) 


no 


r detection 


S-P 


no e production 


/ 


no t production 


strong 


constraints 


chiral supp. 


no 


r detection 


S + P 


no e production 


/ 


no r production 


strong 


constraints 


chiral supp. 


no 


r detection 


T 


no e production 


no P-odd part 


no r production 


strong 


constraints 


chiral supp. 


no 


t detection 



Propagation (/ = e, u, d) 



V -A 


/ 


V + A 


/ 



Table I: Classification of the vertices from Eqs. Q - © according to their impact on reactor and superbeam experiments. 
Terms marked with / can give a sizeable contribution; for all other terms, the reason for their suppression is given (see text 
for details). 



non-pcrturbativc nature of the strong interactions, 
these cannot be easily related to the couplings of 
the fundamental quarks. 

If (V — A)(V + A) couplings to electrons exist, the 
processes in which they appear will in general have 
an energy dependence different from that of the cor- 
responding standard processes. For anti-neutrino 
production in nuclear reactors, however, this dif- 
ference is completely negligible because the neu- 
trino spectrum from nuclear (5 decay is governed 
by kinematical effects and by the Fermi function, 
which describes final state Coulomb interactions. 

The cross section for the inverse (3 decay pro- 
cess, which is used to detect reactor anti-neutrinos, 
is derived from empirical values for the effective 
vector and axial- vector couplings, so any possible 
(V — A) (V + A) contribution is automatically taken 
into account properly. 

Finally, (V — A)(V + A) interactions involving elec- 
trons in the beam detector, are midly chirally sup- 
pressed, in analogy to (V — A)(V + A) interactions 
involving muons. 

10. As we have seen in Eqs. flU) - ©, non-standard 
matter effects can only have a (V — A)(V — A) 



or (V — A)(V + A) Lorcntz structure, as long as 
we restrict the discussion to left-handed neutrinos. 
For the computation of the coherent forward scat- 
tering amplitude, the factor [/7 P (1 ± 7 5 )/] has to 
be averaged over the neutrino trajectory, and for 
unpolarized matter at rest, the only contribution 
is Nf = /7 /, the fermion density appearing in 
the matter potential. Since Nf is independent of 
the axial current, we conclude that both possible 
Lorentz structures would have the same impact on 
the non-standard matter effects. 



To conclude this discussion, we would like to emphasize 
again that non-(V — A)(V — A) Lorentz structures can 
play an important role in reactor and superbeam exper- 
iments. However, these experiments do not have the ca- 
pability to distinguish different Lorentz structures, un- 
less the spectral distortion caused by (V — A)(V + A) 
terms in the superbeam detector is taken into account. 
In the following, we will neglect this spectral distortion 
for simplicity, assuming that it is anyway hidden by the 
systematical uncertainties in the neutrino cross sections. 
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2.3. Hamiltonian approach to non-standard 
interactions in neutrino oscillations 



The obsevrations from the previous section can be ex- 
ploited to further reduce the number of free parameters 
in our problem. To this end, we define effective cou- 
plings e s a p, e d p, and e™p, corresponding to non-standard 
interactions in the production, detection, and propaga- 
tion processes. e s aB describes a non-standard admixture 
of flavor p to the neutrino state which is produced in as- 
sociation with a charged lepton of flavor a. This means, 
that the neutrino source does not produce a pure flavor 
neutrino eigenstate \v a ), but rather a state 



E 

j3= e,/^,r 



'a/3 



(12) 



Similarly, the detector is sensitive not to the normal weak 
cigenstates, but to the combination 



y 0\ 



(M 



E < 

a— e,fj,,r 



a(3\ 



(13) 



Note that in s s aj3 , the first index corresponds to the flavor 
of the charged lepton, and the second to that of the neu- 
trino, while in s^g, the order is reversed. We have chosen 
this convention to be consistent with the literature. 

In general, the matrices (1 + e s ) and (1 + e d ) are non- 
unitary, i.e. the source and detection states are not re- 
quired to form complete orthonormal sets of basis vectors 
in the Hilbert space: 



E 

a— e,/i,r 



E 

j3=e,fjj,T 



4){4\*i, (M) 

* ( 15 ) 



We can read off from Tab. [J that the 3x3 coupling 
matrix e s a p receives contributions from e s ^' d ' V:izA anc [ 

e^p' d ' S±P , while and are built up only from 
(V - A)(V ± A) contributions. 

Since the coefficients e s ea and e d e (for a = e, /i, r) both 
originate from e s ' u ' d ' v±A , the (V - A)(V ± A) coupling 
to up and down quarks, we have the constraint 



e s = e d * 



(16) 



which again reduces the number of independent param- 
eters by 3. The aforementioned spectral distortion and 



mild chiral suppression in the supcrbcam detector could 
invalidate Eq. (|16|) . but we will neglect it in the following. 

Similarly, the (V — A) (V ± A) part of e^ a and e d are 
the same, and since the (S+P)(S±P) Lorentz structures 
have less impact in e 



d aiJL than in e^ a , we will typically have 



\e s I > \e d I 



(17) 



(barring fine-tuned cancellation effects). If we assume all 
non-standard interactions to be of the (V — A)(V — A) 
type, as is sometimes done in the literature, the con- 
straints from Eqs. (|TB|) and (fT7|) are tightened to e s = 

(e d y. 

It is clear from Tab. U that coupling to r leptons is 
irrelevant in our case, so we can also take 



'T/3 







(18) 



for all a, (3. 

e m is an additive contribution to the Mikheyev- 
Smirnov-Wolfenstein (MSW) potential in the flavor basis, 
Vmsw = a cc diag(l, 0, 0), which now becomes 



V 



MSW 



acc 



(\ 














p m 






MM 


°M r 


\ 




fc MT- 





(19) 



with a C c = 2V2G F N e E. Recall from Eq. (TJ} that the 
diagonal entries in this matrix have to be real, so that 
the Hamiltonian will remain hcrmitian, and can be diag- 
onalized by a unitary mixing matrix. 

Since we are interested in a combined analysis of reac- 
tor and superbeam experiments, it is important to keep 
in mind that the effective e matrices are the same for 
both types of experiments, because, under the assump- 
tions and approximations discussed above, those entries 
which may be relevant in both of them and e^ e ) are 
identical in both cases. 

The oscillation probability is obtained as 



=H\e 



d, ~- iHL w s j\ 2 



\(l + e d U (e- lHL )Al + e s ) aS 



7<5 ^ 

'l+e d f e- lHL (l + e s f 



8a 



(20) 



where 



fa 



H a 8 = -rr; 



2E 



u a 



'21 



(21) 



G 



The Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix U is parameterized as 

C12C13 S12C13 s 13 e~ lSc 

U = I -S12C23 - c 12 s 13 s 23 e lScF c 12 c 23 - s 12 s 13 s 23 e l5cr c 13 s 23 (22) 

S12S23 - c 12 s 13 c 23 e l5cr -c 12 s 23 - s 12 s 13 c 23 e l5cF c 13 c 23 



As usual, Sij and Cy denote the sine and cosine of the 
mixing angle 9ij, and Sep is the (Dirac) CP phase. 

For anti-neutrinos, we have to replace e s , e d , and e m 
by their complex conjugates in the above equations, and 
reverse the signs of acc an d Sqp ■ 



2.4. Perturbative calculation of oscillation 
probabilities 

In practice, it is very convenient to expand the oscil- 
lation probabilities in a perturbative series with respect 
to the small quantities #13, Am^/ Am^, and (s^g 1 ' | in- 
stead of attempting to evaluate Eq. ([2"U|) exactly. Fol- 
lowing a procedure similar to the one explained in the 
appendix of [6(| , the first order expansion reads 



+acc 



(0 



U 



(i)t 



A? 



( f! 

'er 



-efx 

.777 

'-/,(> 

fIT 



L 31j 
£ 



(27) 



By C/( 0) , we denote the PMNS matrix for 8 13 = 0, and 
contains the first order terms in #13. The unper- 
turbed Hamiltonian, , can be easily diagonalized ex- 
actly, so that the matrix exponentials in the above equa- 
tions can be evaluated. It is straightforward to extend 
the expansion to higher orders. 



Pu 3 —>vi 



where 



P 



(0) 



P 



(i) 



(23) 



3. MODIFIED NEUTRINO OSCILLATION 
PROBABILITIES FOR REACTOR AND 
SUPERBEAM EXPERIMENTS 



,(0) 













(3a 



(24) 



— iff' ' L 



0a 



0a 

l + e d ) 



e -iHi*L (l + £ sf 



0a 



d\T e _ lH (0) L 





* 






{3a 





dx e 



J (3a 

-ifll »(I-i)|j(l) e -iJJ (I,) a 



h.c. 



0a 

(25) 



and 



JT<°) - A 
^ " 2E 







U<f>) 


':■ ) 




\ A ™3l/ 



[/(0)t 



+acc 



(26) 



^ — 



1 





k \ 


u"» 


Ami 







C/(o)t 



C/(o)t 



Am§ x 



To study the impact of non-standard interactions on 
reactor and superbcam experiments, we need in partic- 
ular the oscillation probabilities Pns^nd, P v s^ rl/ d^ and 

e e jj. e 

-P„s^„d, to which these experiments arc sensitive. We 
will first present approximate analytic formulas for these 
quantities in Sees. 13.11 to 13.31 and then discuss them in 
Sec. 13.41 All approximations were carried out with the 
perturbative method described in the previous section. 
We have checked that the expressions presented in this 
section reduce to the well-known standard oscillation re- 
sults if NSI are absent by comparing them to the expres- 
sions derived in [H, H3, . Moreover, we have verified 
all formulas numerically, term by term, using Mathemat- 
ics. 

To simplify the notation, let us make the abbrevia- 
tions = sin#y, Cij = cos(9,j, s 2X ij = sin 20^-, and 
c 2 xij = cos20jj. Moreover, it will be convenient to split 
the e parameters into their real and imaginary parts by 



s,m,d 

writing e^p 



~'a0 



\exp(i(j> s ^p l ' d ). To keep our results 



as general as possible, we will not impose the constraint 
from Eq. (|16[) . but treat e s and e d as completely indepen- 
dent matrices. Thus, our formulas will be also applicable 
to experiments with fundamentally different production 
and detection processes, e.g. to a neutrino factory, where 
the production occurs through a purely leptonic vve\i 
vertex, while the detection process vlud involves cou- 
pling to quarks. For reactors and superbeams, it is, of 
course, straightforward to impose Eq. (|16[) a posteriori. 
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3.1. The v e 



channel 



In a reactor experiment, the ratio L/E is chosen close 
to the first atmospheric oscillation maximum, so we 
can safely neglect terms proportional to Am^/Amlj. 
Moreover, matter effects arc irrelevant, i.e. we can take 
fl cc = 0. Finally, we will neglect terms suppressed by 
The last approximation implies that 



'13' 



'13' 



we only consider the interference terms between stan- 
dard and non-standard contributions, but not the pure, 
incoherent, NSI effect. This is justified for most real- 
istic extensions of the Standard Model, where e <C 1, 
but it has been pointed out in |28j that, from current 
model- independent experimental limits, the NSI might 
even dominate over the standard oscillations in some sit- 
uations. We find for the oscillation probability 



o i / r , s , . 2 Ami L 

£ eM |si3S23 cos(d C p - (j) ep ) sin 



9 Ami L 
1 -4 S ? 3 sin 2 



+ 2 
-4 

+ 2 

-4 

+ 2 

-4 

- 2 
-4 

- 2 
+ 



£e M l s i3S23 sin(<5 C p 

£erl s 13C23 COs((5 CP 
£erl s 13C23 Sm(<5 C P " 



AE 

1S . AmlL 



sm 



l 31 J 



AE 
Ami L 
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It is interesting to remark that, due to the e s ee and ef e 
terms, this expression can be different from unity even 
for Am^L/AE -C 1, i.e. at the near detector (ND) site. 
Indeed, we obtain in this case 
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where we have taken into account also second order terms 
in e, which may be important in the near detector due 



to the large event rates. Eq. (|2"9"|) corresponds to an over- 
all rescaling of the neutrino flux, which, however, will 
be hard to detect in a realistic experiment due to the 
systematical flux uncertainty. 



3.2. The u u 



v £ channel 



In the derivation of P v s^ v d^ we will relax our approxi- 
mations from the previous section, and take into account 
also terms of 0{s vi Am 2 21 / Am 2 ^), ©([Am^/Am^] 2 ), 
and O (s Am^ / Am^), to reproduce the correct Sep de- 
pendence. For experiments with a relatively short base- 
line, such as T2K, it is justified to assume vacuum os- 
cillations, if the e m parameters are < O(0. 1) (in Sec. 
we will discuss cases where this is not true, and we will 
see that non-standard matter effects can then be large in 
T2K). The vacuum oscillation probability reads 
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The corresponding expression for the near detector is 
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If the baseline is longer, as is the case e.g. in IMOi^A, mat- 
ter effects are important. To keep the notation concise 
in this case, we define the effective 13-mixing angle in 
matter, which is given to lowest order by 
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The oscillation probability is then 
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Most of the ©(Am^/Am 2 ,! ) terms contain factors of Am^/oocj which can be large at low matter densities, and 
might therefore seem to spoil the accuracy of the expansion in this case. However, the oscillatory terms in square 
brackets become small as Am^/acc becomes large, so that overall, the 0(Arri2i/ Am^) terms remain subdominant 
even if the vaccum limit is approached. 



3.3. The v u — > u u channel 



For P^s^^d, we obtain 



pvac _ 1 _ 2 - 2 ^'"31 J 

r vf ^ v d — 1 s 2x23 bm 



in vacuum, and 



AE 

2|e^|cos^ M + 2|e^|cos^ 

Ww\ cos <^ + 2 l4l «»02x23 «n 2 



Am L 

2 ( Kr I COS <^ r + | | COS C 2 x 23 «2 x 23 sill 2 



(l^rl sin ^ r + |e^| sin 52x23 sin A ^ lZ/ 



°(tij;) +0 "'-'> +0 < e2 >- < 34 » 



1/^ — >y u z X Zo a 771 



2|^|cos^ + 2k^|cos^ 



10 



2|^|cos^] S L 2 3Bin 2 Am3lL 



AE 



[2\e 

Am 2 i 

( Kt I sin ^ T + | sin s 2x 23 sin ■ 



|€ T/ J COS T/ J C 2 x23«2x23 sm 



. 2 Am^L 



4£ 



-le m I 
2 i mm' 

5^1 



s 2x23 1 



■ sm ■ 



^ r 2E ~~" 2E 

,2 



2£ 
4s 



2x23C 2x23 co SV ^ S in 



3 2x23 c 2x23 



a C ci . Ato^L 2 a C c . 2 Am§ 1 i 



■ sm ■ 



5 2x23 c 2x23" 



2E 2E 
accL . Am^L 



4 «2x23 c 2x23 



Am§ x 



2£ 



7 -31 

2E 



x 2 a CC . 2 

4s 2 x23C2x23-T— 2-sm 

^"'31 



4# 

Am^L 



AE 



(35) 



in matter. 



3.4. Interpretation 

To understand the physical origin of the formulas de- 
rived in the previous sections, let us consider Figs. Q] -El 
where we show schematically the possible reaction chains 
that a neutrino can follow before its detection. We also 
indicate the respective suppression factors of the transi- 
tion amplitude, but to simplify the discussion and to im- 
prove the clarity of the figures, we do not explicitly show 
contributions proportional to Am^/ Am 2 ^, which would 
appear in concurrence to the (^-suppressed processes if 
013 is very small. Dotted lines indicate suppression due 
to standard effects, while dashed lines represent transi- 
tions that are suppressed by the non-standard parame- 
ters. The thick blue (black) paths are those followed by 
the standard oscillation channels, while light gray lines 
indicate paths that are suppressed by more than one e 
parameter, and are therefore mostly negligible. In Tab.HTl 
we summarize the same considerations in tabular form. 
In the first part of the discussion, we will assume e s and 
e d to be completely independent in order to keep the 
discussion as general as possible. The constraints from 
Eqs. (fT6|) and (fTT]) will only be implemented afterwards. 

On the reactor side, we can read off from Fig. Q] and 
Tab. HH that, in the presence of just one type of non- 
standard interaction, only e s ee , e* , e s erl ef e , e d e , and 
e d e are relevant. Of these, e s ee and e d e have an 0(e) 
effect which could even exceed standard oscillations if it 
were not for the near detector, where these parameters 
would induce a similar effect as in the far detector. Since 
the absolute reactor neutrino flux is not known precisely, 
the measurement relies on the relative counting rates in 
both detectors, so that the impact of e s ee and e d e is can- 
celed. The remaining NSI, £* , e s eT , e d e , and e d e , con- 
tribute to the oscillation probability at 0(£ sin 26*13) (or 
at C(e sin 26*13 + £ 2 ), if the constraint e s ea = e d * e from 
Eq. (fT6|) is implemented). This can be comparable to the 
standard term, so that these NSI are expected to have 



a large impact on the far detector. They do not affect 
the near detector as long as only one of them is present, 
but if Eq. (fTo| is taken into account, the near detector 
will receive an 0(e 2 ) contribution that can be important 
in some situations. All these considerations are nicely 
confirmed by Eq. (|28|) . 

For a superbeam experiment, Fig. [5] and Tab. [TT] show 
that, as long as only one type of NSI is taken into account, 
e^ e , e^ M , e^ r can affect the production process, while 

<4» e^i) <^L> £ U> £ re; and 4^ ma y bc important in 
the detector. The propagation can be affected by all 
entries of e m . As for the reactor case, the suppression 
factors associated with these different processes can be 
understood from simple physical arguments, which are 
confirmed in a more rigorous way by Eqs. (|30p and (|34|) . 
Let us discuss the different types of NSI in more detail: 

• £g M and e d affect only the disappearance channel 
and are therefore irrelevant for the measurement of 
013 (they may, however, lead to wrong results for 
the leading atmospheric parameters). 

• £^ M also affects the disappearance channel, but it 
can also lead to a modified v ^ rate in the near de- 
tector, and therefore to wrong assumptions on the 
initial neutrino flux. This, in turn, could lead to 
a misinterpretation of the far detector appearance 
measurement, so that the $13 measurement is in- 
fluenced indirectly. The effect in the near detector 
is suppressed by e, so it will affect the far detector 
analysis only at the subleading level of £sin 2 26*13. 

• £* M , £* r and e d e are relevant for the appearance 
channel, but the corresponding amplitude is sup- 
pressed by £8^26*13, so that in the oscillation prob- 
ability, we would obtain a subdominant contribu- 
tion of C(e sin 2 26'i3) from the interference of the 
standard and non-standard terms. 

• £ fie- £< fiei anal £ re have amplitudes of 0(e), i.e. 
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e sin 20i3 
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e 
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e sin 20i3 
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"The factor cos 2023 cannot be derived from Fig. [3] but only from Eq. (135 b - 

Table II: Classification of the reparameterized non-standard interactions e s , e d , and e m according to their impact 
transition amplitudes for reactor and superbeam experiments. For each NSI coupling, only the leading order effect is 
The framed entries highlight those terms that are most relevant to the determination of 0i3 (see text for details). 



on the 
shown. 



they contribute to the appearance probability on 
the level of e sin 29\z, which can be comparable to 
the leading contribution ~ sin 2 26 £^ e , however, 
is suppressed by a factor of cos(Am| 1 L/4i?), which 
is small at the first atmospheric maximum around 
which the beam is centered. Note also that the 
modified v e flux in the near detector that is ex- 
pected in the presence of e* e or ef ie can help to 
actually detect the NSI, although part of it may 
be misinterpreted systematical error on the in- 
trinsic beam background. 

• Of the non-standard matter effects, only and 



£™ contribute at leading order to the appear- 
ance probability -P™^^. Of these, e™ is already 

strongly constrained experimentally (47| . and so is 
not expected to have a large impact on reactor and 
superbeam experiments, e™, on the other hand, 
could contribute significantly to the superbeam ap- 
pearance channel, in accordance with [1, [3(3 ■ All 
other non-standard matter effects are suppressed 
by an additional power of S13, (or, more correctly, 
S13, which is, however, still small since we are far 
from the MSW resonance). 

It is interesting to observe that the sensitivity to 
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Figure 1: Possible contributions of e s and e d to the event rate 
in the far detector (a), and the near detector (b) of a reactor 
Ve disappearance experiment. Thick blue lines indicate the 
reaction chain for standard oscillations. Dotted lines indicate 
processes that are suppressed by standard three-flavor effects 
proportional to 613 or Amii/Amjj!, while dashed lines rep- 
resent transitions that are suppressed due to non-standard 
interactions. Paths which would only be accessible in the 
presence of two different non-standard effects, are shown in 
light gray since they are usually subdominant. 




(c) 



is very weak, although this type of interaction 
corresponds to a simple rescaling of the standard 
MSW potential. However, it is not a leading order 
effect, and therefore does not appear in our approx- 
imate formula, Eq. (133[) . 

• In p™ at the dominant matter effect is e"i, and 

since there is no #13 suppression from the inter- 
ference with the standard amplitude, this effect is 
even stronger than those in p™ at Note that, 

from Fig. [31 one might expect £^ and £™ to be 
of similar strength as e™ r , but when one performs 
the calculation, it turns out that an additional sup- 
pression factor C2x23 appears (cf. Eq. ([55]) ). 

• The implementation of the constraints e s ea = e^* 
(Eq. (H) and |e*J > |e*J (Eq. does not 
lead to any new effects, except for the appearance of 
an additional e 2 term in the disappearance channel 

for e s = £ d * 

Let us finally emphasize the crucial importance of the 
standard and non-standard phases in the oscillation prob- 
abilities: The formulas from Sees. 13. 11 - 13. 3l reveal that un- 



Figure 2: Possible contributions of e s and e d to the event rate 
in a superbeam experiment for the appearance channel (a), 
the disappearance channel (b), and in the near detector (c). 
The meaning of the colors and line styles is the same as in 

Fig.m 

favorable phase combinations may suppress non-standard 
effects, even if the modulus of the corresponding e param- 
eter is large. 



4. SIMULATION OF REACTOR AND 
SUPERBEAM EXPERIMENTS 

To fully assess the high-level consequences of non- 
standard interactions for realistic reactor and superbeam 
experiments, we have performed numerical simulations 
using the GLoBES software [69l [70l] . We have considered 
the following scenarios: 

• T2K + Double Chooz. Our simulation of the T2K 
far detector, Super-Kamiokande, is based on [7l| . 
Most parameters are taken from the T2K letter of 
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(a) 




(b) 



Figure 3: Possible contributions of e m to the event rate in the 
superbeam appearance channel (a), and in the corresponding 
disappearance channel (b). The meaning of the colors and 
line styles is the same as in Fig. [1] 



intent [48|, an d the systematical uncertainties are 
based on [72|]. We include a separate 1.0 kt wa- 
ter Cerenkov near detector with similar properties 
as the far detector, and similar systematical un- 
certainties. To model the interplay of the two de- 
tectors, we introduce a common 10% uncertainty 
on the neutrino flux, and a common 20% error 
on the number of background events in the v e ap- 
pearance channel. In the absence of non-standard 
interactions, these correlated errors would cancel 
completely, since the total neutrino flux and the 
background contribution are effectively calibrated 
by the near detector, but if s s ' d ^ 0, this calibra- 
tion can be wrong, and there may be an observable 
effect. The neutrino interaction cross sections are 
taken from [6(| [zH ■ We assume 3 years of neutrino 
running and 3 years of anti-ncutrino running, each 
with a beam power of 0.77 MW. The fiducial far de- 
tector mass is 22.5 kt, and the baseline is 295 km. 
We consider v e appearance events as well as the 
disappearance signal. The background for the dis- 
appearance channel is made up of neutral current 
events, while for the appearance measurement, neu- 
tral current events, misidcntified muons, and the 
intrinsic beam backgrounds can contribute. 

For the simulation of Double Chooz, we use the 
same parameters as in [531 ] . and the cross sections 
for inverse beta decay are taken from [23]. As for 
T2K, we simulate the near and far detectors sepa- 
rately, but take into account the appropriate cor- 
relations between systematical errors. In particu- 
lar, we introduce a 2.8% flux normalization error, 
which is correlated between the near and far de- 



tectors, uncorrelated 0.6% fiducial mass errors for 
both detectors, uncorrelated 0.5% energy calibra- 
tion uncertainties, and an 0.5% bin-to-bin uncorre- 
lated error. 

• NO^A + DC-200, where DC-200 refers to a reac- 
tor experiment similar to Double Chooz, but with 
a 200 t far detector [53|. Such a large reactor ex- 
periment has a considerable sensitivity not only to 
the total event rate, but also to distortions of the 
energy spectrum. 

The simulation of the v e appearance signal in NOz^A 
is based on [E(j, while for the disappearance 
channel, we follow [2H]. We assume 3 years of neu- 
trino running and 3 years of anti-neutrino running, 
with a beam power of 1.12 MW. The far detector 
mass is 25 kt, and the baseline is 812 km, with 
an average matter density of 2.8 g/cm 3 along the 
trajectory, while the near detector has a mass of 
0.0204 kt, and is located at 1 km from the target. 
Again, we introduce, in addition to the uncorre- 
lated systematical errors from [5(J [75| , a correlated 
10% uncertainty on the total neutrino flux, and a 
correlated 20% error on the v e background. 

The parameters and systematical errors of the DC- 
200 scenario are identical to those of Double Chooz. 

Unless indicated otherwise, we calculate the respective 
event rates using the following "true" values for the os- 
cillation parameters [76| : 

sin 2 20f 2 ue = 0.84, 



sin 2 26>^ uo 
sin 2 26»f 3 ue 



1.0, 
0.05, 



£cp c = 0.0, 



(36) 



(Am 2 ,! 
(Am 2 



CP 
2 \truc 



31 



7.9 x 10" 5 cV 2 , 
| truc = 2.6 x 10" 3 eV 2 , 



and assume a normal mass hierarchy. To analyze the 
simulated data, we follow the statistical procedure de- 
scribed in the appendix of [7l[, and define the following 
X 2 function 1 



channel bin r , j , , 

X =mm ^ ^ 



true -true 



)-N ij (X,e = 0)\' 



iV ij (A true ,e truo ) 

+ Priors, (37) 

where Nij denotes the number of events in the i-th 
energy bin for oscillation channel j, the vector A = 



1 In the implementation of superbeam experiments, we assume 
the events to follow the Poisson distribution. However, for il- 
lustrative purposes, it is sufficient to consider the more compact 
approximative Gaussian expression. 
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(#12; $i3i #23) <5cPj Am^, Am| 1; 6) contains the standard 
oscillation parameters and the systematical biases b, and 
e represents the non-standard parameters. In the fit, 
we marginalize x 2 over all standard oscillation parame- 
ters and over the systematical biases, but since we want 
to study how a standard-oscillation fit gets modified if 
there are non-standard interactions, we keep the NSI pa- 
rameters fixed at 0. 2 The prior terms implement ex- 
ternal input from other experiments and have the form 
(x — x true ) 2 /<7 2 , where x stands for any oscillation pa- 
rameter or systematical bias, and <j x is the correspond- 
ing externally given uncertainty. We assume 6\2 to be 
known to within 10%, and A m 2 ,-, to within 5% from so- 
lar and reactor experiments |76| . When analyzing the 
reactor experiment alone, we additionally assume a 15% 
uncertainty on #23 and a 5% error on Amf^. Beam ex- 
periments are themselves sensitive to #23 and Am§ 1; so 
we omit these priors for them. 



5. NSI-INDUCED OFFSETS AND 
DISCREPANCIES IN 13 FITS 

Using the simulation techniques discussed in the previ- 
ous section, we can now determine the errors that are in- 
troduced when non-standard interactions are present in 
reactor and superbeam experiments, but are not prop- 
erly taken into account in the respective fits. Possible 
outcomes of such fits are shown in Fig. H] for our two 
scenarios. As "true" parameter values, we have taken 
sin 2 2#i3 = 0.05 and <5cp = 7r, and the NSI contribution 
was assumed to be e™ = 0.5e~ l7r / 2 in the upper panel, 
and e s eT = e d Te = 0.05 (fulfilling Eq. (HI) in the lower 
panel. These NSI parameters are rather large, but still 
consistent with current bounds [3l|, According to 

the discussion in Sec. [3J only the superbeam experiment 
should be affected in the first case, while in the second 
case, there should be an impact on both experiments. 
The shaded areas show the 90% confidence regions for 
the reactor experiment, while the colored contours are 
for the superbeam. The data has been calculated under 
the assumption of a normal mass hierarchy, but the fit 
has been performed for both the normal mass ordering 
(solid blue/dark gray contours) and for the inverted or- 
dering (dashed pink/light gray contours). The vertical 
black line and the colored diamonds represent the re- 
spective best fit values, while the black star stands for 
the assumed "true" parameter values. We have assumed 
two degrees of freedom for the superbeam experiments, 
and one degree of freedom for the reactor setups, which 
are insensitive to Sep- 

It is obvious from the plots that the standard oscil- 
lation fit to 9 13 and <5cp can be severely wrong if NSI 



2 Of course, when computing confidence intervals for certain pa- 
rameters, we have to keep these parameters fixed as well. 



are present. This observation is similar to the one made 
in [34| for a neutrino factory. In the case shown in the 
upper plots of Fig. 01 the reactor experiment gives the 
correct best fit value, but the superbeam results conflict 
with this measurement. In the case of NOvA, we even 
obtain a fit value above the Chooz bound and a fake hint 
to the mass hierarchy, indicated by the fact that, within 
the resolution of the plot, the 90% contour reduces to a 
single point. In the second case (lower plots), both ex- 
periment agree very well, but they erroneously seem to 
rule out the "true" #13. 

Of course, the NSI scenarios analyzed in Fig. [H were 
only two examples, and a more systematic analysis of 
non-standard interactions in reactor and superbeam ex- 
periments is desirable. This is done in Figs. O - [71 
where we show how the (standard oscillation) #13 fits 
in T2K/Double Chooz respectively in NOi^A/DC-200 may 
be distorted in the presence of non-standard interactions. 
For each diagram, only one of the independent e param- 
eters was assumed to be non-zero, but we have ensured 
that Eqs. (fl^|) and (TTT|) are fulfilled. In particular, we did 
not consider the hypothetical case e* a = 0, e£ ^ 0, but 
onl y ^ a + 0, e d ali = and e* Q = ^ 0. Moreover we 
omit all entries of e m except , because they are either 
strongly constrained already (e^), or do not have any 
impact on the #13 measurements (e™, , e™, and e™ T ). 
The modulus of each parameter has been varied between 
and its current upper bound, which is 0.1 for e s ^ from 
universality in charged lepton decays [3l|, and 0.7 for 
[47J . 3 The complex phases were allowed to vary be- 
tween and 2ir. For each combination of |e| and arg(e), 
we have then performed a fit assuming standard oscilla- 
tions, and the resulting best fit values for #13 are shown 
in the plots. Points giving a good fit (better than 3er in 
both experiments), are drawn as thick colored lines, with 
the hue indicating the respective value of |e|. Dark red 
(dark gray) corresponds to |e| = 0, while yellow (medium 
gray) corresponds to the upper bound of |e|. Points giv- 
ing a fit quality worse than 3(7 arc shown by thin light 
gray lines, and the information on |e| is omitted for them. 
All computations have been performed for two different 
"true" values, sin 2 26>i 3 = 0.01 and sin 2 29 13 = 0.05, as 
indicated by the black stars. 

By comparing the plots with Tab. [Til we find that our 
expectations for the impact of the different s parameters 
from the discussion in Sec. [3J arc confirmed. A partic- 
ularly interesting situation arises for ef e , because this 
parameter has a sizeable effect in both, the reactor ex- 
periment and the superbeam setup. It is especially dan- 
gerous because it induces a similar offset in both exper- 
iments, i.e. one would find perfectly consistent #13 fits, 



3 Note that, according to the naive estimate from Eq. J7ji such 
large values of |e| would correspond to Mnsi ~ O(100 GeV). In 
many models, such low new physics scales are already ruled out, 
but in our model-independent treatment, they are still viable. 
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T2K / Double Chooz (90% C.L.) NOvA / DC-200 (90% C.L.) 




0.01 0.02 0.05 0.1 0.01 0.02 0.05 0.1 

Sin 2 26> 13 Sin 2 26>i 3 



Figure 4: Two examples for the errors that are introduced if non-standard interactions are neglected when fitting #13 and Sep 
to the data of reactor and superbeam experiment. In the upper plots, a discrepancy arises between the two experiments (the 
NOi^A fit is even above the Chooz bound), while in the lower plots, there is a common offset, leading to consistent results, but 
erroneously "ruling out" the true #13 (indicated by the black star) at a high confidence level. The left hand plots are for T2K 
and Double Chooz, and the right hand ones are for NO^A and DC-200. The gray shading represents the 90% confidence region 
from the reactor experiment, and the vertical black line shows the corresponding best fit value for #13. The 90% contours from 
the superbeam are shown as solid blue (dark gray) lines for a normal hierarchy fit, and as dashed pink (light gray) lines for 
an inverted hierarchy fit. The colored diamonds represent the corresponding best fit values. In interpreting the computed \ 2 
values, we have assumed 2 degrees of freedom for the beam experiments, and one degree of freedom for the reactor setups. 



which might, however, be far away from the true value. 

Other parameters may lead to fit points far from the 
diagonal, corresponding to seemingly conflicting fits. An 
interesting case is the e^ e term, for which the non- 
standard interaction mimics a significantly modified v e 
flux in the near detector. This, in turn, leads to a mis- 
calibration of the beam- intrinsic backgrounds, so that, 
at the far site, many of the actually oscillation-induced 
v e events will be mistaken as background. Thus, the fit 
value for #13 becomes too small. However, we can also 
read off from the plot, that, in this situation, the quality 



of the standard oscillation fit becomes so bad that the 
NSI effect can actually be detected. Note that the curves 
for large |e^ e | look slightly untidy, because for some pa- 
rameter values, the smallest x 2 is provided by the normal 
hierarchy fit, while for others the inverted hierarchy fit is 
marginally better. Therefore, frequent "jumps" between 
these two solutions occur. 

When interpreting Figs. [5] and [5J it is important to 
keep in mind that the error bars of the experiments con- 
sidered here are rather large (cf. Fig. [4|), so that even 
sizeable deviations from the diagonals will in most cases 
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Figure 5: Distortion of the #13 fits in T2K and Double Chooz in the presence of e s a p and £% a . For each plot, the modulus of the 
corresponding e parameter has been varied from to 0.1, and its phase from to 2ir. For every such combination, we show the 
result of standard oscillation fits of #13. Connected lines represent contours of equal |e a a| and varying phase. Dark red (dark 
gray) lines correspond to \s a &\ = 0, while yellow (medium gray) lines correspond to \e a p\ = 0.1. Points giving a quality of fit 
worse than 3<r in at least one of the two experiments are plotted in light gray. The black stars indicate the assumed "true" 
sin 2 2#i3. 



only create some tension, but no unambiguous contradic- 
tion between the beam and reactor fits. 

Of the non-standard matter effects, we expect from 
Tab.HTlthat only e™ and e™ should have any effect on the 
#13 fits, e™ is already strongly constrained from charged 
lepton flavor violation experiments [47j , but e™ may still 
have a large impact. In fact, for extreme values of this 
parameter, there is even the possibility that NOi'A would 
erroneously report a #13 value above the Chooz bound. 

Let us emphasize that, in order to obtain reliable es- 
timates for the impact of non-standard interactions on 
reactor and superbeam experiments, it is crucial to take 
the information from the near detectors into account. To 



show this, we have also studied how Figs. [5] and [6] get 
modified if we use a simplified simulation, in which the 
near detector does not appear explicitly, but only through 
suitably small values for the systematical uncertainties. 
In doing so, we have again treated e s and e d as com- 
pletely independent matrices, possibly violating Eqs. (JTHJ) 
and (TPT|) . Thus, the results are also applicable to setups 
where e s and e d are indeed unrelated. In accordance with 
our expectations from Tab. [Til we have found: 

• For e s ee and ef e , the effect on the reactor becomes 
stronger without the proper treatment of the near 
detector, because these terms no longer cancel then. 
Moreover, the discovery potential becomes worse, 
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Figure 6: Effect of ^ and £p a on the 0i3 fits in NO^A and DC-200. The color-coding is the same as in Fig. [5] 



i.e. there will be fewer gray segments in the plot. 

• For the superbeam, the discovery of e^ e and £™ e 
becomes also much harder without the near detec- 
tor, because the clear signature of an apparently 
modified v e flux at the near site is no longer avail- 
able. Moreover, for e^ e , the strong impact of the 
NSI on the #13 fit in the superbeam, which we have 
identified as a near detector effect in the above dis- 
cussion, vanishes in the single-detector simulation. 

• The contours for e s eT = e^* are deformed without 
the near detector because in this (unrealistic) sit- 
uation, it is no longer possible to misinterpret the 
non-standard effect as a reactor flux calibration er- 
ror. Note that this misinterpretation is only due 

are identical, since the 



to the fact that e s eT and e 



near detector is only affected if both are present 
(cf. Tab. HI)) , Otherwise, it would retain its capa- 
bility to properly calibrate the reactor flux to its 
true value. 

So far, we have only considered situations in which one 
non-standard parameter is dominant, and all others are 
negligible. In realistic models, however, many parame- 
ters may be of the same order of magnitude. Since it 
is impossible to visualize the resulting high-dimensional 
parameter space, we resort to the scatter plots shown 
in Fig. [HI These plots were created by choosing a ran- 
dom value for each NSI parameter, and then performing 
a standard oscillation fit to the resulting experimental 
data. The moduli of the e parameters were logarithmi- 
cally distributed between 10 -8 and their current upper 
limits, where we have assumed the model-independent 
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Figure 7: Effect of e™ on the 613 fits in T2K and Double 
Chooz, resp. in NOf A and DC-200. We do not show plots for 
the other entries of e m , since these are either strongly con- 
strained already (e^), or do not affect the #13 measurement 
( £ eej e ™T, and £??-)■ The color-coding is the same as in 
Fig. [5] but the scale is different since the bound on is 
weaker than that for e^'» [3ll, |47|. 



bound from universality in charged lepton decays [3l| for 
e s and e d , and the results of [47] for e m . The phases were 
distributed linearly between and 27r. 

We can see from Fig. [5] that there are again points 
which yield a clear discrepancy in the #13 fits of the reac- 
tor and supcrbcam data, and others which correspond to 
a common offset of the fit value. The color coding shows 
that for a considerable fraction of the parameter space, 
the non-standard effect can actually be discovered. It 
is interesting to observe that there are some points for 
which the reactor fit lies above the Chooz bound. This 



indicates that already with the present data, some parts 
of the parameter space could be ruled out. 



6. DISCOVERY REACH FOR NON-STANDARD 
INTERACTIONS IN A COMBINED ANALYSIS 
OF REACTOR AND SUPERBEAM DATA 

Let us now discuss the prospects of actually detecting 
the presence of non-standard interactions in reactor and 
superbeam experiments. We define the discovery reach as 
the range of e parameters for which the quality of a stan- 
dard oscillation fit is below a given confidence level. In 
Figs.[9]to[T2l we show numerical results for this quantity, 
which were obtained by performing standard oscillation 
fits to the combined data of both experiments. 

The results can again be interpreted with the help of 
Tab. HIl and of the formulas derived in Sec. [31 We see that 
for those non-standard parameters which have a large im- 
pact on any of the observed oscillation channels, there is 
typically also a good discovery potential. For some pa- 
rameters, it comes from the reactor measurement, for 
others, it is dominated by the superbeam. It is remark- 
able, however, that in the case e s eT = ef*, there is prac- 
tically no discovery potential at all, because neither ex- 
periment can discover these parameters on its own, and 
there is also no significant discrepancy between them, but 
only a common offset in their #13 fits. 

It is interesting to observe that the good discovery 
reach for e* r comes from the disappearance channel, as 
can be easily verified from the corresponding analyti- 
cal formulas in Sec. [3j Note that in those plots where 
= e^* is assumed, there is no discovery potential 
because the corresponding NSI terms in Eq. (jM)) cancel. 

The discovery reach depends strongly on the phases 
of the NSI coupling constants, 0*'n' m . To first order in 
S13, all off-diagonal entries of the e matrices (except the 
and e^ T contributions in the — * disappearance 

channel) are accompanied by a combination of </>^'™ and 
5<jp. To first order in Am^/Am^, they typically appear 
together with factors of cos </>^g or sin 0^ . The diagonal 
components of the e matrices usually have prefactors of 
cos 02?^. 

The plots do not exhibit any phase dependence in the 
discovery reach for e* e and e^ e because the sensitivity 
to these parameters comes mainly from the modified v e 
flux in the near detector of the superbeam experiment. 
We have checked, that, in accordance with Eq. ([50)) . the 
phase dependence would reappear if the near detector 
were omitted in the simulation. 

Turning to non-standard matter effects described by 
e m , it is clear that the discovery potential will be very 
limited, since already standard matter effects are small 
in T2K and NOz^A, and completely negligible in Double 
Chooz and DC-200. Therefore, we use a different scale 
for the horizontal axis in Figs. Ql] and [T2l However, for 
some entries of e m , the present bounds are very weak. 
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Figure 8: Possible outcomes of standard three- flavor oscilla- 
tion fits to reactor and superbeam data in the presence of non- 
standard interactions. Each plot contains two datasets, one 
for sin 2 2fl t 1 r 3 uc = 0.01, and one for sin 2 26\ r 3 uo = 0.05 (indicated 
by the black stars). For each dataset, 5,000 random combina- 
tions of s parameters where chosen, with their moduli being 
distributed logarithmically between 10~ 8 and the respective 
upper bounds [3lL H3] , and their phases varying linearly be- 
tween and 2-7T. Each of these random non-standard scenarios 
was then fitted under the assumption of standard three-flavor 
oscillations. Green (dark gray) points indicate a very good 
quality of this fit, while light gray points denote a fit qual- 
ity worse than 3<r in at least one of the two experiments, i.e. 
an effective discovery of the non-standard effect. The plots 
show that non-standard interactions can induce ostensible dis- 
crepancies between reactor and superbeam data (off-diagonal 
points), or a common offset (close-to-diagonal points), which 
would lead to consistent, but wrong results. Note that some 
points lie even above the Chooz bound. A reactor fit above 
the Chooz bound indicates that the corresponding combina- 
tion of NSI parameters and #x3 Ue could already be ruled out 



and s™, but they cannot compete with the current 



In particular we have < 1.0, < 0.7, and < 
1.4 [ij. Figs. [TT] and [fi] thus show that the bound on 
could be improved by NO^A, but not by T2K. We 
should, however, keep in mind that, according to Eq. ((TJ), 
|e™| ~ 0.7 corresponds to M NS i ~ 100 GeV, and it is 
hard to imagine a model, that could yield such a low 
NSI scale without violating present electroweak precision 
data. Both T2K and NO^A have some sensitivity also to 

e; 

bounds < 5 ■ 10~ 4 and e™ < 0.1 

Note that, according to Eqs. ([53"]) and (f3l>|) . the sen- 
sitivity to e™ and e™ comes from the v e appearance 
channel, while the sensitivity to £™ r has its origin in the 
disappearance channel. 

Let us dwell for a moment on the interesting shape of 
the sensitivity contours for e™ and e™ , which can only 
be understood by taking into account terms proportional 
to |e| 2 . Let us consider, for example, e™. According to 
Eq. (|33|) . the NSI contribution to the oscillation proba- 
bility is, to first order in |e| and neglecting Am^, 
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By carefully studying this expression, one finds that the 
energy dependence in the first and third terms of Eq. (j3"5)) 
is quite different from that of standard oscillations, which 
is proportional to sin (Am 3I — acc)L/AE. Therefore, 
these terms will be easy to detect, while the second term, 
which modulates the spectrum in the same way as stan- 
dard oscillations, can be absorbed into a modified #13 fit, 
and will therefore be hard to detect. From the phase de- 
pendence of these terms, we expect that, for our choice 
of S^p c = 0, the discovery reach should be good for 
(j)™ ~ 0, 7r, and poor for 0™ ~ g7r, §7r. The plots in 
Figs. [Tl] and [12] reveal that the discovery reach indeed 
shows this behavior, except for an unexpectedly good 
sensitivity at <j>™ = |7t. To understand this, we have to 
take into account the second order terms, which we have 
found to be 
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The important observation is that the net effect of the 
second order terms is always positive, while for the first 
order terms, it is positive at <f>™ ~ |-7r, and negative at 
</>™ ~ \-k. In the first case, we would therefore need a 
much stronger deviation of the fitted #13 from its true 
value in order to absorb the non-standard term. This, 
however, is disfavored by the reactor measurement, so 
that the combined fit improves the discovery reach by a 
considerable amount at 0™ — §vr. We are here in the 
interesting situation that the combination of seemingly 
redundant data sets can be beneficial if there are devi- 
ations from standard three-flavor oscillations. For most 
other non-standard parameters, the discovery reach is 
dominated by either the reactor or the superbeam. 

7. CONCLUSIONS 

In this paper, we have studied the impact of non- 
standard neutrino interactions on upcoming reactor and 
accelerator neutrino experiments. We have first classi- 
fied the allowed NSI terms in the Lagrangian according 
to their Lorentz structure, and have found that many 
of them are irrelevant to reactor and superbeam se- 
tups. Those which can have an impact are mostly of 
the (V — A)(V ± A) type, but in superbeam experiments, 
also (S + P)(S± P) type effects can be important. Since 
reactor and superbeam experiments are not able to dis- 
tinguish different Lorentz structures, we have reparam- 
eterized the NSI coupling constants in order to greatly 
reduce the number of free parameters in the problem. 

Using this reparametcrization, we have then derived 
approximate analytic expressions for the non-standard 
neutrino oscillation probabilities, both in vacuum and 
in matter of constant density. We have developed an 
intuitive understanding of the terms relevant to specific 
oscillation channels, and have classified them accordingly. 

In the second part of our work, we have performed de- 
tailed numerical simulations using GLoBES. We have con- 
sidered two scenarios: T2K combined with Double Chooz, 
and NOz/A combined with a 200 t reactor experiment, 
dubbed DC-200. Our simulations take into account pa- 
rameter correlations, degeneracies, and systematical er- 
rors, and in particular, we employ a realistic treatment 
of the near detectors. We have found that non-standard 
interactions can have a sizeable impact on future reactor 
and superbeam experiments, if the coupling constants are 



close to their current upper limits, and if complex phases 
do not conspire to cancel them. The biggest impact is 
on the 9\3 measurement: If NSI are not properly taken 
into account in the fit, the results may be significantly 
wrong. There are scenarios in which a clear discrepancy 
between reactor and superbeam experiments shows up, 
but we can also have the situation that both fits sets seem 
to agree very well, but the derived #13 value has a signif- 
icant offset from the true value. It is even possible that 
the true #13 is erroneously "ruled out" at 3er. To detect 
this kind of problems, a third experiment, complemen- 
tary to the other two, would be required. Thus, we see 
that the possibility of non-standard effects should always 
be kept in mind when planning or analyzing upcoming 
experiments. 

We have also studied the discovery potential for NSI 
in reactor and superbeam experiments, i.e. the range of 
non-standard parameters, which can actually be detected 
by these experiments because the quality of a standard 
oscillation fit becomes poor. We have found that, de- 
pending on the complex phases, some NSI may be dis- 
covered if their coupling constants are not more than a 
factor of 5 smaller than the current upper bounds. The 
best discovery reach is obtained only if both, reactor and 
superbeam experiments, and also the respective near de- 
tectors are considered in the analysis. In most cases, 
one of the experimental channels dominates the discov- 
ery reach, but there are also situations where only the 
discrepancy between the single-experiment fits indicates 
the presence of NSI. Our discussion thus shows that re- 
actor and superbeam measurements, which might seem 
to be redundant in the standard three-flavor framework, 
turn out to be highly complementary once non-standard 
effects arc considered. 
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